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Abstract. We consider systems of diffusion equations that have considerable interest in 
■ Soil Science and Mathematical Biology and focus upon the problem of finding those forms 

of this class that can be linearized. In particular we use the equivalence transformations of 
the second generation potential system to derive forms of this system that can be linearized. 
In turn, these transformations lead to nonlocal mappings that linearize the original system. 
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> ! 1 Introduction 

. Whilst systems of pure diffusion equations, in both their linear and nonlinear forms are well 

known and have many physical and biological applications, the research described here focusses 
' on less familiar cases where diffusion coefficients or other 'shape' functions are defined either in 

■ general or poor analytic terms. Of particular interest here is the case of the extension of Richard's 

equation, which describes the movement of water in a homogeneous unsaturated soil, to cases 
describing the combined transport of water vapour and heat under a combination of gradients 
Ch ■ of soil temperature and volumetric water content. Such coupled transport is of considerable 

^ . significance in semi-arid environments where moisture transport often occurs essentially in the 

^> I water vapour phase [S]. Under these conditions the transport equations, valid in a vertical 

^ • column of soil, may be written in the form [S] 
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(1) 

where u{x,t) and v{x,t) are respectively the soil temperature and volumetric water content at 
depth X and time t. It is important to realize that extensions which include the coupled diffusion 
of solute follow in an obvious way. 

If we introduce the potential variable w, we can generate the auxiliary system of ([T|). 

= u, 

wt = f{u, v)ux + g{u, v)v^, 

vt = [h{u,v)ux + k{u,v)voi:]^ . (2) 
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Introduce potential variables w and z we generate the second generation auxiliary system of , 
= u, 

wt = f{u, v)ux + g{u, v)vx, 

zt = h{u,v)ux + k{u,v)vx. (3) 

The study of Lie symmetries of system has been considered in |21 E]- It is known that 
Lie symmetries of the auxiliary systems, in some cases, lead to non-local symmetries known as 
potential symmetries, for the original system. Lie symmetries for the systems ((2| and (jHI that 
induce potential symmetries for the original system have been considered in |^. 

Here we calculate the equivalence transformations for the systems Q-®- We use the equi- 
valence transformations of Q to derive those forms of (jSJ that can be linearized. Consequently 
these point transformations lead to contact transformations that linearize the corresponding 
forms of Furthermore for one case of 0, which admits infinite-dimensional Lie symmetries, 
we construct a linearizing mapping. 



2 Equivalence transformations 

An equivalence transformation is a nondegenerate change of the independent and dependent 
variables taking a PDE into another PDE of the same form. For example, an equivalence 
transformation of the system will transform it into a system of the same form, but in 
general, with different functions f{u,v), g{u,v), h{u,v) and k{u,v). The set of all equivalence 
transformations forms an equivalence group. We use the infinitesimal method [7j to determine 
the equivalence transformations of the systems (|T|)~-(j3|). 

2.1 Equivalence transformations for system 

We seek for equivalence group generators of the form 

t^F ^ d ^ d d d ,5 , d . d ,9 
ox at ou ov Of og oh Ok 

We will present the results, without giving any detailed calculations. The method for determining 
equivalence transformations is presented in [2]. (See also in j^). We find that the system has 
a continuous group of equivalence transformations generated by the following 10 infinitesimal 
operators: 

yE _ 9 vE _d vE _ d vE _ 9 

d 9,9 d ,9 9 ,9 

= - ^ "« • = ^ "8/ <* - ^' 8^ - 

d 9 / „ , N 9 9 
^- = ^9^-^97 + ^^-'%+^9a:- 

Using Lie's theorem we show that the above equivalence transformations in finite form read 



x' = ClX + C2, t' = C^t + C4, u' = C5U + CqV + C7, v' = CgU + CqV + ClQ, 
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where 

cj [05(09/ - cgg) + CQ{cgh - csk)] , _ cf [-C5(c6/ - c^g) - ceiceh - c^k)] 



C3(C5C9 - CeCs) ' C3(C5C9 - CeCs) 



^/ ^ cf [Cgjcgf - Csg) + Cgjcgh - Csk)] ^ ^ Cj [-CfjjcQf - c^g) - Cgjceh - Csfc)] 

C3(C5C9 - CeCg) ' C3(C5C9 - CgCs) 

Clearly, if the functions /, g, h, k are linearly dependent, then one of the functions of the 
transformed equation can be taken equal to zero, provided that the constants involved satisfy 
certain relations. That is, system can be mapped into a system of the same form but with 
one of the functions to be zero. For example, \i f = g — h + k then the mapping x ^ x, t ^ t, 
u ^ u + V, V ^ u + cv, c ^ 1 transforms the system into a system of the same form, but 
with (7 = 0. 

2.2 Equivalence transformations for system ((SJ) 

We find that the system @ has a continuous group of equivalence transformations generated 
by the following 9 infinitesimal operators: 

d d d d d „d d ,5 ,9 

ou ov ow ox ow oj og on ok 

^8 =*di~^df~%~^dh~ dk' ^9 -^a^-^^s^ + '^s^- 

The equivalence transformations, in finite form, read 

x' = c\x + C2, = est + C4, n' = C5U + cg, = C7t; + cs, 

= CiCsW + CiCgX + Cg, 

where 

/' = ElZ 5' = ^5cfg ^, ^ C7cf/t _c\k 
ca ' C7C3 ' C5C3 ' C3 

2.3 Equivalence transformations for system 

We find that the system Q has a continuous group of equivalence transformations generated 
by the following 14 infinitesimal operators: 

7E _d yE _d yE _ d yE _ 9 

d d d 

+ (-ufc + 37;g + n/) — + vh— + (2t;fc + uh) — , 
Of? an, ok 

+ + {vk + 3n/i -vf )^ + {2uk - vg) 

d d d d ,(9 ,9 

ax aii aw a/ og On ok 
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Naturally using Lie's theorem we can find the equivalence transformations in finite form. We 
find that 

x' = ax + piw + P2Z + ci, t' = —t + 5, w' = qix + aiw + a2Z + C2, 

7 

z' = q2X + a^w + 042 + C3. 

However the forms of u' and v' cannot be found easily from Lie's theorem. In the next section 
we show that they take the form 

, qi + aiu + a2V , q2 + a^^u + 04?; 
u = , V = . 

a + piu + p2V a + piu + p2V 

Furthermore we state that the forms of g' , h' and k\ which are all linear in /, g, h and k, 
are very lengthy. 

In the following section we use the finite form of the equivalence transformations to derive 
special cases of the system that can be linearized. In turn these transformations lead to 
nonlocal mappings that linearize the corresponding forms 

3 On linearization 

Here we consider the problem of finding forms of the nonlinear system ^ that can be linearized. 
We adopt the idea of the transformation 

x' = v, t' = t, u' = -, v' = x (4) 

u 

that maps the auxiliary system of the nonlinear diffusion equation ut = [u~'^Ux]x 0) 

Vx = U, Vt = U~'^Ux 

into the auxiliary system of the linear diffusion equation u'^., = u'^,^,, 

Vx'=u', v[,=u'^,. 
The above transformation is a member of the equivalence transformations of 

Vx =U, Vt = f{u)Ux. 

Motivated by the above results we find the forms of /(n, v), g{u, v), h{u, v) and k{u, v) such that 
system Q can be mapped into a linear system by the equivalence transformations admitted 
by ©. These local mappings will lead to nonlocal mappings that linearize the corresponding 
forms of the system 

We consider system (Q) with the four functions equal to constants. That is, it takes the linear 
form 

u't' = fj'iu'x'x' + f^2Vx>x'^ = Hsu'x'x' + ^^'^Vx'x'■ (5) 
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We introduce the potential variables w' and z' to obtain the auxiliary system of © 



u 



4' = /^3<' +/^44'- 



(6) 



From the equivalence transformations of the system (jS)) we deduce the finite transformations 

x' = ax + piw + P2Z , w' = qix + aiw + a2Z, z' = q2X + asw + a^z (7) 

where we have taken, without loss of generality, the translation constants equal to zero. Clearly, 
the inverse transformations is of the form 



X = hx + riw' + r2z' , w = six' + biw' + b2z' , z = S2x' + h^w' + h^z' . 
Also from the equivalence transformations we deduce that 
t = -ft', 

with the translation constant taken to be zero. 

Introducing vector notation, we can write systems Q and © in the form 



(8) 



and 



u, wt = F{u)u^ 



u 



w. 



respectively, where 
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u = 
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F(u) 



f{u,v) g{u,v) 
h{u,v) k{u,v) 



Furthermore transformation ^ can be written in the form 
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and the inverse transformation (IHl) takes the form 



(9) 



(10) 
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^1 






. ^3 


/i4 _ 



(11) 
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where 



Therefore from the transformations and 1)12(1 we deduce that 
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(12) 



(13) 
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where I2 is the 2x2 identity matrix. In addition: 



b " 




V T 1 

a p 




ab + r ■ q 






" 1 " 


s B 




.1 A . 




sa + Bq 


sp^ + BA 




h 



(14) 



Using normal transformation rules: 
dx' dt' 



Wn 



OX ox 



and from the transformations and (|12jl becomes 

Wx = {s + Bw'^,) {a + p- w^) . 
Thus from the two potential systems © and (jlOj) 

It = (s + Bu') {a + p ■ u) . 

Finally, from the relations and (|14() we obtain 
Au + (7 



Hence, 



u 



a + p ■ u 



qi + aiu + 



q2 + a3« + 04^; 



These forms of u' and u' can also be obtained from the equivalence transformations, however, 
in a more complicated manner. 
Consider now 



I dx I dt 



and from the transformations and p2p and the two potential systems Q and (|1U|) we find 
Aw^, = (q + ^M)(r^A<.,) + ^AFux. (15) 
Multiply by r'^ and using the relations and (|14|1 it follows that 
-iP^Fux 



a + p ■ u 



and also consider 

, . dx , dt 

which gives 

u'^i = (6 + r^u')u'^ 
and using the above form of u' and the relations (|14|1 reduces to 



It, 



a + p • u 
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These expressions simplify ()15|) to 

Ku'^ = -^{q + Au){p^ Fux) + -i{a + p ■ u)AFux. 

Using the relations (|T1^ and (fTl)) and differentiate the expression of u' with respect to x the 
above equation takes the form 

A [(a + p^u)Aux — (p • u^)Au — pu^q] = —7(9 + Au){p'^ Fu^) + 7(0 + p ■ u)AFux 

and simplifying to get 

A [(a + p^u)A - {Au + q)p^] = 7(0 + p^uf [{a + p^u)A - {Au + q)p^] F. 

Solving for F{u) to obtain 

[H{u)Y'^KH{u) 



F{u) 



(16) 



7(a+p^it)2 ' 
where 

H{u) = (a + p^u)A - {Au + q)p^ . 
Summarizing we have: 

Theorem 1. The nonlinear system of diffusion equations Q can he mapped into the linear 
system ()10() by the equivalence transformation admitted by @ if and only if the functions F{u) 
is of the form (|16)) . 

Remark 1. Point transformations that linearize system (j^j) lead to contact transformations 
that linearize system (^. 

Remark 2. In the case where A = I2, that is, system (jSJ becomes two separate linear diffusion 
equations with diffusivity constants equal to 1, the linearizing form of is 



ut 



{P1U + P2V + a)2 



vt 



(17) 



{P1U+P2V + a)2_ 

It is known that the part of transformation ((3}, namely x' = v, t' = t, v' = x, which is known 
as pure hodograph transformation maps the potential equation vt = v~'^Vxx into the linear heat 
equation v'^, = v'^,^,. In the spirit of the work in where such transformations were classified 
for the potential equation vt = f{vx)vxx, we consider the potential system of 

Wt = f{Wx, Zx)Wxx + g{Wx 1 Zx ) Zxx ) 

zt = h{wx, Zx)wxx + k{w 

XI Zx)Zxx- 

(18) 

Transformations that presented in this section which linearize systems of the form ^ , can also 
be employed to linearize systems of the form (|18|) . We present the results in the following 
theorem: 

Theorem 2. The nonlinear system of potential diffusion equations (|18|) can be mapped into the 
linear system 



fJ-lWx'x' + lJ'2Zx 



fJ-S'^x'x' + f^4:Zx'x' 



by the transformation ^ if and only if the functions f{wx,Zx), g{wx,Zx), h{wx,Zx), k{wx,Zx) 
are of the form 



f{wx,Zx) g{w xi Zx) 

h{Wx,Zx) k{Wx,Zx) 



-f{a + p^Wxy 



where 



H{wx) = (a + p^Wx)A - [Awx + q)p^ 
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Using the above theorem the corresponding form of Remark 2 reads: 
Remark 3. The nonlinear potential system 



wt = 7 : : — vTi 



{piWx + P2Zx + a)^' {piWx + P2Zx + ay 

can be transformed into two separable linear heat equations. 

4 Examples of linearizing mappings 

In this section we use the results of the previous section to present two examples. 

Example 1. We consider the system (jlTf) with its corresponding second generating potential 
system 

Ux 

Wx = U, Wt 



{piu + p2V + ay ' 

Zx = v, zt = — (19) 

{P1U + P2V + ay 

We have shown that system (|19|) can be linearized using the equivalence transformations of ©. 
That is, the transformation 

x' = ax + piw + P2Z, t' = t, w' = qix + aiui + a2Z, z' = q2X + asui + 042;, 

, qi + aiu + a2V . q2 + a^u + a^v 

u = , V = 

a + piu + p2V a + piu + p2V 

maps the linear system 

w'x' = u', w[, = 4, = v', z[, = 

into the nonlinear system (|19|1 . 

This transformation leads to the contact transformation 

dx' = adx + (n + v)dx + \{piu + p2V + a)~'^{ux + Vx)\ dt, dt' = dt, 

, qi + aiu + a2V , q2 + asu + a^v 

u = , V = 

a+piu + p2V a + piu + p2V 

that maps the two separate linear diffusion equations 

into the nonlinear system (|17)) . Using Theorem [21 we deduce that the potential form of (|17|) . 

'U^xx Zxx 
Wt = 7 ■ ■ TTT, Zt 



{piWx+P2Zx + ay' {piWx + P2Zx + ay 

can be linearized by the transformation 

x' = ax -\- piw + P2Z, t' = t, w' = qix + aiw + a2Z, z' = q2X + asw + a^z. 

Example 2. We consider the special case of the equivalence transformation of @ 

, , , u+ 1 , , 

X =w + z, t =2t, u = , V = , w = x + z, z =x + w 

u + v n + -y 
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which maps the hnear system 
into the nonUnear system 



Wx = u, Wt 



(n + v)^ 



Ux + 



-1^1 + (1^2 - 1^3)'" + 1^4 
(u + 



1^1 "W^ + (1^2 — J^sl'W — 1^4 —l^lUV + 1^3^ + 1^2''^ — 2^4 
=7;, = ; : Ux H -, Vrr 



{u + vY 



where 



1^1 = m + tJ-2 - fJ-a - M4- 

Clearly if we set = ^4 = 1, i-t2 = fJ'S = 0, we obtain a special case of the previous example. 
Now the above transformation lead to the contact transformation 

dx' = (n + v)dx + [i'i{vux — uvx) + i'2{ux + Vx)] {u + v)~^dt, dt' = 2dt, 

, 1 + f , 1 + u 
u = , V = 

U + V U + V 

that maps the linear system 



Uf = fJ-lUx'x' + f^2Vx'x', 

into the nonlinear system 

ViUV + 1^2 + I'sv + 1/4 



= fJ'Su'x'x' + ^^'iVx'x' 



ut 



vt 



{u + vY 



-UlV? + {V2 - V^)U + 1/4 
'^x H 7 i 



Ux H : 



5 A linearizing case of the system ([2]) 

In Section 3 we used the equivalence transformations of © to derive linearizing mappings. We 
point out that employment of the equivalence transformations of Q does not lead to linea- 
rizing mappings. However such linearizing mapping, which is not member of the equivalence 
transformations, exists for a special case of @. 
We consider the special case of © 



Wx = U, Wt = -U "^Ux, Vt=[u ^Vx], 

which is the first generation potential system of 



(20) 



Ut 



[U V],, Vt= [U ^Vx]^. 



(21) 



System (|2Up admits the infinite-dimensional Lie symmetries 
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where the function (j)(t, w) satisfies the backward Unear heat equation (f>t + (j)ww = and ^l^{t, w) 
satisfies the hnear heat equation ipt — ipww = 0. These symmetries induce potential symmetries 
for the system ()21() . 

If a nonlinear PDE (or a system of PDEs) admits infinite-parameter groups, then it can be 
transformed into a linear PDE (or into a linear system of PDEs) if these groups satisfy certain 
criteria. These criteria and the method for finding the linearizing mapping using the infinite- 
dimensional symmetries can be found in 1^. Hence, using the method described in [3], the above 
infinite-dimensional Lie symmetries of (|2nj) lead to the transformation 



1 

X' = w, t = t, w = v = V, w 

u 



t' = t, - - ' 



which maps the linear system 



into the nonlinear system (|2U|1. Consequently this mapping lead to the contact transformation 

/ 1 



dx' = tidx -|- u ti^di, di' = dt 
which maps the linear separable system 



^t' — ~'^x'x' ' ^t' — '^x'x' 



" n = — , V = V 



into the nonlinear system H21|). 

We point out that the above result cannot be achieved using the equivalence transformations 
of © . However it is straight forward to derive the above contact transformation as a special case 
of the one derived in the Example 1, Section 4, which maps the two linear separable diffusion 
equations u'^, = u'^,^, and v'^, = v'^,^, into the nonlinear system ut = and vt = ['>J'~'^Vx]^- 



6 Conclusion 

In this paper we have considered the problem of finding forms of the general class of systems of 
diffusion equations Q that can be linearized. We have employed the second generation potential 
system © and derived the equivalence transformations admitted by this system. Using these 
transformations we classified special cases of (jSJ that can be linearized. In turn these transfor- 
mations lead to nonlocal mappings that linearize the corresponding forms of Furthermore 
we determined a special case of that can be linearized, by considering the Lie symmetries of 
the first generation potential system ((2)), which induce potential symmetries for the system (|T]). 

The question that arises here is: Are these the only cases of the system Q that can be linea- 
rized? Furthermore: Can the results obtained here be generalized for systems of n equations? 
These are two of the problems of our investigation in the near future. 
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